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Yukawa effects on the clock onboard a drag-free satellite 
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ABSTRACT 

The Yukawa correction to the Newtonian gravitational force is accepted as a parame- 
terization of deviations from the inverse-square law of gravity which might be caused 
by new physics beyond the standard model of particles and the general theory of rela- 
tivity. We investigate these effects on the clock onboard a drag-free satellite: dynamics 
of the satellite and influence on the time transfer link. We find the Yukawa signal 
in the time transfer is much harder to detect with current state of clocks than those 
effects on the dynamics, especially the secular change of periastron, by laser ranging 
in the case of an artificial Earth satellite carrying a frequency standard with an orbit 
of a = 10 7 m and e = 0.01. 
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1 INTRODUCTION 

It is gravitation that is the first of the four known fundamen- 
tal forces to be understood quantitatively and empirically. 
However, unlike the strong, weak and electromagnetic inter- 
actions, it is also gravitation that still can not be explained 
as a result of the quantum exchange of virtual bosons. The- 
oretically, despite some candidate theories proposed in re- 
cent years, unification of gravitation and the rest part of 
physics remains as a grand challenge. Experimentally and 
observationally, searching for poss ible deviations fr o m the 
inverse-square law of gravity [see lAdelberger et all (|2003l ) 
for a review] , which are predicted by some possible schemes 
of quantum gravity, might shed light on new physics. One 
of the parameterizations of the deviations is the Newtonian 

ravitational potential with a n additional Yukawa correction 

Fischbach et alJll986l Tl992h that is 



V = V N (r) + VykW, 

where 

t r / \ Gm 1 m,2 

V N (r) = , 

r 

Gmim 2 
VYK(r) = aexp 



(1) 

(2) 
(3) 



Here a is a di mensionless strength paramete r and A is a 
length scale [see iFischbach fc Talmadgel (|l999h for a review 
of constraints on a and A] . 
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Inspired by the idea of orbital tests of relativistic grav - 
ity using Earth satellites l|Damour fc Esposito-Fareselll994h , 
we would like to investigate the Yukawa effects on the clock 
onbo ard a satellite and their observabi lity. In this methodol- 
ogy (|Damour fc Esposito-Farese 1 19941') . the sa tellite will use 
drag-free technology (jPughll 19591 ; lLangdll964h . making it be 
able to follow a free fall trajectory, and its three-dimensional 
motion is tracked at the centimeter or even less level either 
by laser ranging or the global positioning system. The ad- 
vantage of using an space-borne clock is that Yukawa ef- 
fects, in principle, will appear not only in the dynamics of 
the satellite but also in the time transfer signals. 

In Section [2j treating the Yukawa correction as small 
disturbance, we will follow the standard procedure of per- 
turbation in the celestial mechanics and obtain the Lagrange 
plan etary equations of t he sa tellite. It is worth mentioning 
that lHaranas fc Ragosl (|201ll ) derived the same equations 
but their results contradict the conservation law of the an- 
gular momentum in the central force problem (see Sec. [2] 
for details). We will study the Yukawa effects on the time 
transfer of the onboard clock in Section [3] In Section |4j the 
observability of these effects is discussed. Conclusions and 
discussion will be debriefed in Section 



2 YUKAWA EFFECTS ON SATELLITE 
DYNAMICS 

There is no doubt that an artificial Earth satellite prac- 
tically feels a lot of non-gravitational disturbance, such as 
atmospheric drags and radiation pressures. These effects are 
usually larger or much larger than the relativistic corrections 
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to the Newtonian gravity and some effects caused by possi- 
ble deviations from the inverse-square law of gravity. A very 
promising way t o suppress them to a ti ny level is the drag- 
free technology (Pugh 1959; Langc 1964) that is already em - 
ployed by the Gravity Probe B (GP-B) jEveritt et al.H201ll ) 
and Gravity fiel d and steady-state Oce an Circulation Ex- 
plorer (GOCE) (|Drinkwater et al.ll2003l ). Its basic principle 
is that a space vehicle offsets the nongravitational orbital 
disturbances by staying around an inside test mass that is 
shielded from the environment. Since the test mass is gov- 
erned by the gravitational force only, the satellite is able to 
follow a free fall trajectory in the background gravitational 
field. In the following parts of this work, we assume that the 
clock used to detect the Yukawa effects will be carried on- 
board such a drag-free satellite so that we focus on the pure 
gravitational interactions. Since the relativistic celestial me- 



dBrumbere & 


Kooeikinl 1 19891: ISoffell 


198S 


D amour et al. 


1994 iKoDeikin et alj 


2011 



), we will pay at- 



tention to the Yukawa effects only. 

With widely used notations in celestial mechanics, if a, 
e, i, £1, uj and M are taken as independent orbital arguments, 
the e volution of a satellite orbit under a perturbing potential 
R is (|Danbvlll962h 
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In the case of Yukawa correction to the Newtonian potential, 
according to Eq. l[3]l. 



R = — a exp 
r 



(10) 



where = GM. The above expression shows R depends on 
r only, i.e. R = R(r), which means the two-body problem 
under the framework of the Newtonian gravitation with the 
Yukawa correction is a special case of the central force prob- 
lem. Evaluating the partial derivatives involving R, we have 
the Lagrange planetary equations for the satellite as 



dM 
~dT 



(1 



nae r 

2 u 

-\ 7 -aexp 

na 2 r 



e 2 ) n 

j a exp 



(16) 



It is worth mention ing that, in the work by 
lHaranas fc Ragosl (|201ll ). the authors derived the same 
equations but their results are quite different from above 
ones. Especially, they found the orbital inclination i varies 
due to the existence of t he Yu kawa force [Eqs.(18) and 
(33) in lHaranas fc Ragosi l|201ll )]. Given the fact that the 
Newtonian gravity with the Yukawa correction is a special 
case of central force so that the angular momentum per 
unit m ass of the 2-body system h is conserv ed as a constant 
vector (|Danbvlll962l ; iGoldstein et al.ll2002D . the inclination 
i = arccos(h. • z/|h.|), where z is the unit vector along z-axis 
of the inertial reference system, will never change. 

For secular evolution of the orbit, we need to average R 
over one orbital revolution of the satellite that is 



R 



an 

Rdt — — exp 
a 



(17) 



where Ip(z) is the modi fied Bessel function of the first kind 
( Arfkcn fc Weberj[2005h . Therefore, we can have the long 
term variations of the orbital elements as 
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(18) 
(19) 
(20) 
(21) 
(22) 

(23) 



where h{z) = dl {z)/dz |Arfken fc Weberll2005h . Our re- 
sults Eqs. II17P and (|22[) identically match those given by 
lloriol (|2012l ). With the dynamics of the satellite under the 
Yukawa perturbation, we can study its effects on the time 
transfer of the onboard clock. 
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(13) 
(14) 



l + -)cos/, (15) 



3 YUKAWA EFFECTS ON TIME TRANSFER 

For a satellite surrounding the Earth, the Geocentric Coor- 
dinate Time (TCG) is the coordinate time to describe its 
dynamics and the time variable in its equations of motion 
dSoffel et al.ll2003l ; iKopeikin et"aill201ll ; ISoffel fc Langhans! 
120131 ). However, the actual reading of a clock onboard the 
satell i te is the proper time r llMisner et al ]| 19731; ISoffel et all 
120031 ; IKopeikin et ail 1201 ll ; ISoffel fc Langhansll2013D . They 
diverges with time. To trace the readings of a space-borne 
clock to a clock on Earth or other time scales, the transfer- 
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mation between TCG and r is usually taken as an interme- 
diate. 

Accordin g to t he relativistic theory of tim e scales 



ilSoffel et al.l 120031; iKovalevskv fc Seidelmarml |2004 



iKopeikin et alj l201ll ; ISoffel fc Langhansl l2013h . the trans- 
formation between r and t = TCG is 

r-B 



At 



1 ,r 

1 + ^ + 2^ 



dr + 0{c 



where, if taking the Yukawa correction into account, 



1 + a exp 



(24) 



(25) 



In the above expression of V, the contributions from the 
Earth's oblateness and external bodies are not metioned and 
we assume they can be properly modelled to adequate pre- 
cision. (We will go back to this issue in the next section 
of obs ervability.) Following the procedures used by I Nelson! 
|201ll ). we can obtain 
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(26) 



in which E is th e eccent r ic ano maly. When a = 0, it returns 
to the result of iNelsonl (|201lh . From the Eq. IpB"]). we can 
see that the Yukawa effects on the time transfer have two 
components: a secular term and a periodic term. The secu- 
lar term might not be observable because it is mixed with 
the intrinsic frequency drift of the clock. The periodic term, 
whose frequency is as same as the leading term due to gen- 
eral relativity, might show up, but its observability depends 
on its amplitude 



-4, 



fiae exp 



1 + 



(27) 



4 OBSERVABILITY OF YUKAWA EFFECTS 

This section will be dedicated to an important issue: the 
observability of these Yukawa effects. We focus on a relative 
small astronomical length scale that the space vehicle stay 
in the vicinity of Earth and take the domain of Yukawa 
parameters as 

V = {(A,a)|10 6 m < A < 10 10 m, 10~ 12 < a < 10~ 8 } (28) 
based o n the literature [su ch as iFischbach fc Talmadgel 



Il999l): IXdelbereer et al.l feood); iDeng et al.l 
Lucchesi fc Peronl (|20l6T )7 iDend (|201lh ; iLucchesi 



2009) 



2011) 



Ioriol ( 2012T )]. We adopt a satellite with a = 10' m and 
e = 0.01 as an example. 

According to Eq. (|22p . the domain T> can generate 
the secular change of uj, (uj), ranging from several nano- 
arcsecond per year to hundreds arcsecond per year (see color 
indexed Fig. [TJ . If the (lj) can be determined to the level of 
micro-arcsecond per year by laser ranging that the current 
best fitted values can reach better than ~0.1 milli-arcsecond 
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Figure 1. Color indexed secular change of u according to Eq. 
([22} in V. 



per year IjLucchesi fc Peronll2010l : lLucchesJ|j)llh , most re- 
gions of the Fig. [1] will be brought into the detectable thresh- 
old. 

For the time transfer, the amplitude Ay [Eq. (|27[) ] of 
the periodic term due to the Yukawa correction holds the 
value from ~ 10 -24 s to ~ 10 -16 s in the domain T> (see 
color indexed Fig. [2}. Even if a current optical clock with a 
fract ional frequency inaccuracy of 8.6 x 10 -18 (|Chou et al.l 
l2010h would be carried onboard the satellite, it is hardly 
able to detect the Yukawa signal in the duration of one or- 
bit. Since picking up a periodic term usually needs to sam- 
ple the data in the time interval comparable with its period, 
which is ~ 10 4 s in this case, such a long time will make the 
uncertainty of the clock's frequency exceed the Yukawa sig- 
nal. One of possible ways to improve the situation with the 
given state of clocks is to dig the higher frequency periodic 
terms (but perhaps with much less amplitudes) if exist. 

At the frequency stability level of a few parts in 10 18 , 
a very important issue involved in the time transfer is the 
contributions [in Eq. I|25|l] from the Earth's oblateness and 
the gravitational tidal fields of external bodies, at least for 
the Moon an d the Sun. They have b e en well inves t igated and 
formulated inlBlanchet et all (|200lh . lSoffel et all (|200St ) and 
iPetit fc Luzuml l|2010l ) and they need to be carefully taken 
care in order to pick up the Yuakawa signals. 



5 CONCLUSIONS AND DISCUSSION 

In this work, we investigate the Yukawa effects on the clock 
onboard a drag- free satellite, including the dynamics of the 
space vehicle and the time transfer of the clock. We find 
the Yukawa signal in the time transfer with current state of 
clocks is much harder to detect tha n those in the dynamics , 
especially the secular change of w l|Lucchesi fc Peronl l2010l : 
iLucchesllioTlh , by laser ranging in the case of a satellite 
with a — 10 7 m and e = 0.01. 

However, some open issues remain. One of them is 
whether the Yukawa effect will reach the detectable realm in 
high frequency terms in time transfer. After all, Eq. (|26|l is 
truncated and other terms do exist but perhaps with much 
less amplitudes. Another interesting issue migh t be these e f- 
fects on the interplanetary time transfer link jDend l2012h. 
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Figure 2. Color indexed amplitude Ay according to Eq. 1271 in 
V. 

which need to step into the region with a larger astronomical 
scale to length. 
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